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Abstract 

We introduce a new method for selecting exclusive semileptonic charmless -B-decays 
in the presence of a large background. The method can be applied to charged and 
neutral 5-mesons decaying into any exclusive neutral or charged hadronic final 
state. The method is designed for high luminosity e + e~ colliders operating at the 
T(4S). It employs an improved partial reconstruction technique for .D*-mesons and 
a novel 0-C event fit to both 5-meson's decay products resulting in the kinematics 
of all particles (including neutrinos) in the event. The charged lepton energies are 
accessible from 1.0 GeV to the kinematic limit. 



1 Introduction 



Semileptonic B decays into charmless final states provide information about 
1 14ft | . Although the relation between \ V u b\ and the exclusive charmless semilep- 
tonic branching ratio is model dependent, this situation is improving as lattice 
calculation results become available [1]. In measuring the decay rates, it is de- 
sirable to sample as much of the kinematically allowed phase space as possible. 
This increases signal efficiency and allows to differentiate between theoretical 
models. 



Branching fractions of exclusive semileptonic B decays into non-charm final 
states have been measured[2] to be approximately 2 x 10~ 4 . Because this small 
rate leads to a relatively large background from b — > c decays and other 
sources, we have developed[3] a new method with excellent background rejec- 
tion and acceptable signal efficiency. The method is based upon the kinematic 
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reconstruction of the entire event with one 5-meson decaying into the low 
rate channel of interest (the "signal") and the other .B-meson into a channel 
with a large branching fraction (the "tag"). The method can only be used at 
e + e~ colliders that operate at the T(AS) center-of-mass energy because the 
kinematic reconstruction needs the magnitude of the £?-meson's momentum 
as an input. If the beams have a non-zero crossing angle or unequal energy 
the method applies after a transformation of all observed final state particles 
to the e + e~ center of mass. The method can be used for charged and neu- 
tral _B-mesons decaying into any exclusive neutral or charged hadronic final 
state. The method requires large luminosities such as will become available 
at S-Factories. The choice of a tag's decay channel is dictated by the need 
to have a large branching fraction and a good detection efficiency The decay 
B — > D*tvi has the largest branching fraction but it has the disadvantage of 
introducing a second neutrino in addition to the neutrino on the signal side. 
As discussed in Sec. 2, events containing two exclusive semileptonic B decays 
can be reconstructed if all particles' 3-momenta, except the two neutrinos, are 
measured. 

To maintain efficiency, the D* is reconstructed using an "improved partial re- 
construction" technique that, like the "standard partial reconstruction" tech- 
nique, only uses the pion kinematics from the decay D* — > Dn, see Sec. 3. 
The method is further improved by an novel 0-C event fitting procedure, in- 
troduced in Sec. 4. 



2 Event Reconstruction 

We consider the reconstruction of events of the type B — > Xi + v^ (signal), 
B — > Yl~V t (tag), where the B, B, X, and Y may be charged or neutral. 
We assume that the momenta of the X, £ + , Y, and l~ are measured. Particle 
masses are assigned according to the signal hypothesis. Incorrect assignments 
will generally cause the event to fail the procedure. Thus the energies of the 
four measured particles are known. The unknowns are the two angles associ- 
ated with the direction of the back-to-back BB and the momenta of the two 
neutrinos, a total of 8. Energy-momentum conservation for each B decay gives 
8 relations so the kinematics of all particles in the event (including the two 
neutrinos) can be calculated. Solutions come in pairs due to square-roots in 
the equations. 

A geometric construction provides further understanding of the method and 
demonstrates its powerful rejection of background events owing to three re- 
quirements, see Fig. 1. The endpoint of the momentum vectors of the back-to- 
back 5-mesons lie on the surface of a sphere of radius 325 MeV/c centered at 
their production point O (the 11 BB sphere"). Starting with the B signal side, 
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Fig. 1. Geometric construction that illustrates the event reconstruction of 
B -» X£+i/ t and B -► WF*. 

the sum of the momenta of X and i + is calculated and shown as the vector 
OP. Energy conservation applied to the B decay gives the neutrino energy 



E Vl = Eb eam — yPx + rn 2 x — yp%+ + mj. The neutrino momentum vector must 

lie on a sphere of radius E v centered on point P. Constraint # 1: require that 
this sphere and the BB-sphere intersect. The intersection is a circle perpen- 
dicular to OP on which the B momentum vector is constrained to lie. On the 
tag side, the sum of the momenta of Y and £~ is calculated and shown as the 
vector OQ. The point Q is the center of a sphere with radius E Ve , calculated 
from energy conservation applied to the B decay. Constraint # 2: require that 
this sphere and the BB-sphere intersect. The intersection is a circle perpendic- 
ular to OQ on which the endpoint of the B momentum vector is constrained 
to lie. Constraint # 3: require that the B and B momentum vectors, whose 
endpoints are constrained to lie on their respective circles, be back-to-back. 
To find the solution, reflect one of the circles through point O and find the 
intersection of the reflected circle and the other circle. If the circles intersect 
at all, there will be two intersections. Background events will generally not 
satisfy all three constraints. 

To maintain acceptable efficiency for signal events, the D* is not reconstructed 
from its Dn decay products: applicable branching fractions would cause an 
event loss of an order of magnitude. Instead the D* kinematics are inferred 
from the kinematics of the it using the "partial reconstruction" technique. To 
have the best possible determination of the D* kinematics (the other three 
particles are measured with negligible errors) for signal efficiency and back- 
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ground rejection we developed an improved partial reconstruction technique, 
to be discussed next. 



3 Partial Reconstruction 



Because the D*-D mass difference is barely large enough to create a pion, 
the decay products in the D* — > Dn decay have small momenta in the D* 
center-of-mass system. This feature is exploited in the long known partial 
reconstruction technique [4] to get a approximate D* momentum using only 
the pion's kinematics. Information on the D is not needed. 

The energy in the laboratory of the D* is given by 

= rnp*E n 1 
D * E* l + (3*(3 D *cos9* y } 



and the angle 9 between the D* and the tt by 



tan/) 2g;sinfl*l+/%fo.cosfl* 



Variables with a * are to be evaluated in the D* center-of-mass system. In the 
standard partial reconstruction technique one sets f3* — 0. This gives 

m D *E n 

Ed* - -jp (.6) 



and 9 = 0. A better approximation is obtained by studying the two-dimensional 
histogram of Ed* and E n , see Fig. 2, for the case B — > D* + £~Ve, D* + — > 
D°7i + . The ISGW[5] model was used to simulate the decays. The dots are 
the average D* energy in each bin of p n . The line represents Eq. (3). There 
is a significant difference between the dots and the line, especially at larger 
energies. 

A lookup table can be constructed to obtain Ed* from E n (or p n ). We show in 
Fig. 3 the difference between the true E D * and E D * obtained with the standard 
partial reconstruction technique (solid histogram) and the improved partial 
reconstruction technique (dashed histogram) that uses such a lookup table. 
The improvement is clear. The improvement can be seen even more clearly 
in the difference between the true pd* and the partially reconstructed Pd*, 
see Fig. 4 The rms deviation of this distribution decreased from 475 MeV/c to 
330MeV/c as a result of the improvement. The distribution of E D * energy in 
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Fig. 2. Two-dimensional histogram of Ed* and E n in B° — > D* + £~V e , D* + — > D 7r+. 
The dots are the average D* energy in each pion momentum bin. The line represents 
Eq. (3). 

a given E v bin is very skewed and for some applications the most probable 
instead of the mean may be more appropriate. A similar method may be used 
to improve the estimated direction of the D*. We found that the benefits from 
improving the D* direction are small. 

In principle there is model dependence in the improved partial reconstruc- 
tion technique. We have compared two models, ISGW[5] and ISGW2[6] that 
both describe the B° — > D* + £^u e and found the difference between them to 
be negligible relative to their respective difference from the standard partial 
reconstruction. 



4 Event Fit 



Because the D* momentum from improved partial reconstruction has an un- 
certainty of approximately 330MeV/c, signal events may fail the event recon- 
struction. An example is shown in Fig. 5 where the solid sphere centered on 
point Q/ is the result of the improved partial reconstruction of the D* while the 
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Fig. 3. Difference between true D* en- 
ergy and the D* energy from partial 
reconstruction. The dashed histogram 
is for the standard and the solid his- 
togram is for the improved partial re- 
construction. 
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Fig. 4. Difference between true D* mo- 
mentum and the D* momentum from 
partial reconstruction. The dashed 
histogram is for the standard and the 
solid histogram is for the improved 
partial reconstruction. 
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Fig. 5. The same event as in Fig. 1 showing that it now fails the event reconstruction 
because of the use of the partial reconstruction technique for the D*. Point Q (Q/) is 
the endpoint of the vector sum of the lepton and the true (partially reconstructed) 
D* momentum. 

true D* leads to the dashed sphere centered on point Q. This is the same event 
as shown in Fig. 1. Now the event does not satisfy constraint ^2 and would 
therefore be rejected as being background. It is possible to recover such signal 
events by varying the D* momentum (causing Q/ and E-y to vary as well) such 
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that an intersection of the sphere centered on Q/ and the BB sphere exists 
(constraint # 2) and that appropriate back-to-back Bs exist (constraint # 3), 
see Sec. 2. This corresponds to performing a 0-C fit. Of all D* momenta that 
satisfy the three contraints, we choose the one where the variation of the D* 
momentum is minimal as measured by a x 2 5 defined as follows. 

Using a model to simulate B — > D*iv , D* — > D tt we measured the covariance 
matrix Vij defined as 

V H = Jj 12 (Ppr,i ~ Ptr,i)(P P r,j ~ Ptrj) h 3 =X,y,Z (4) 

events 

where p tn (p P n) is the i-th component (i = x,y,z) of the true (partially 
reconstructed) D* momentum. Because there is significant dependence of the 
covariance matrix upon the pion momentum, it is measured in bins of pion 
momentum. A x 2 is defined as 

X 2 = 12 12(Pi ~ Ppr^iV'^ijiPj ~ Pprj) hj =X,y,Z (5) 

i 3 

where Pi is the i-th component of the D* momentum. The D* momentum 
with the lowest \ 2 an d satisfying the three constraints is selected. For signal 
events, Xmin i s typically smaller than 1.0 while background events generally 
have very large xLin- 

In Fig. 6 we show the difference between the true and the fitted D* momentum 
(solid histogram). For comparison, we also show the difference of the true D* 
momentum and the D* momentum from respectively the improved (dashed 
histogram) and the standard (dotted histogram) partial reconstruction. Each 
histogram has the same number of events. The fit significantly improves the 
D* kinematics, and that of its associated neutrino, thereby validating the fit 
procedure. 



5 Results 

The event reconstruction technique has been evaluated using simulated events. 
The simulation uses the ISGW model to describe exclusive semi-leptonic B 
decay and the detector simulation of the CLEOII detector [7]. The simulation 
is known to describe the detector response well. The simulation generates raw 
data in the same format as real data. The simulated raw data are processed 
through the same analysis procedures as the real data. For definiteness we 
limit the discussion that follows to the reactions B° — > ix~ t^v^ (signal), B — > 
D* + £-V e (tag) with D* + -> D°ir+. 
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Fig. 6. Difference between the true and the fitted D* energy (solid histogram). Also 
shown are the difference between the true D* momentum and the D* momentum 
from respectively the improved (dashed histogram) and the standard (dotted his- 
togram) partial reconstruction. 



A detailed discussion of the analysis is outside the scope of this paper. Briefly, 
events with at least 4 charged tracks are selected, two of which must be lep- 
tons with momenta greater then 1.0 GeV/c. The leptons can have opposite or 
equal charges because of B°-B mixing. At least one of the other two tracks 
must have a momentum between 40 and 200 MeV/c, the kinematically allowed 
range for the pion from D* + decay. The most important background process 
is the one in which the signal B also decays into a charmed final state. To 
suppress this background at least partly, we calculate the effective mass and 
the total momentum of unused tracks and energy deposits in the electromag- 
netic calorimeter. For signal events all unused tracks and energy deposits are 
the result of D° decay so if their effective mass is greater than the D° mass 
the event is rejected. Likewise, if the angle in the laboratory between their 
total momentum and the pion momentum is greater than 37 deg, the event is 
rejected. 

With this event selection, we loop over all combinations of four tracks in 
an event and over all events and apply the 0-C event fit and the improved 
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partial reconstruction. The resulting Xmm distribution is shown for signal and 
background events in Fig. 7. The latter are BB events that decay according 
to known decay modes, not including semi-leptonic charmless decays. When 
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Fig. 7. Xmin distributions for simulated signal (top) and background (bottom) events. 
Note the logarithmic scale and the number of overflow events in each case. 



we require xLm — 1A on ^J 19 ou ^ °f §53 background events pass, a rejection 
factor of 45. 

The signal efficiency e, defined to include the branching fraction B(D* + — > 
D°7r + ) = 0.68, is measured to be 3.66%. The number of signal events is 
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expected to be 



N = 21 — — — i 3— ^ e 6 

Ufb"V I 1.0 x 10- 4 ) y 1 



where L is the integrated luminosity. The numerical factor includes the B — > 
D* + £~T7g branching fraction and the BB production cross section. 5-factories 
aim for integrated luminosities of at least 30fb _1 /year. A branching fraction 
B(B° — > n~£ + V£) = 2 x 10~ 4 [2] would result in about 50 events per signal 
channel with about 20 events background. When several signal channels are 
combined these numbers are expected to give measurements of \V u b\ whose 
precision is dominated by theoretical uncertainties. A vertex detector would 
improve the tracking of low momentum particles (important for partial recon- 
struction) and allow the rejection of some of the remaining charm background. 
As indicated earlier, the method can be applied to charged and neutral -Bs 
decaying into any exclusive neutral or charged hadronic final states, allowing 
many measurements of \V u b\. 
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